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^: 

^ , We present the general formalism required to derive generalized parton distri- 

ly^ ' butions within a convolution model where the bare nucleon is dressed by its virtual 
meson cloud. In the one-meson approximation the Fock states of the physical nucleon 

(N ; 

, are expanded in a series involving a bare nucleon and two-particle, meson-baryon, 
. 

1^ I states. The baryon is assumed here to be either a nucleon or a A described within 

■ the constituent quark model in terms of three valence quarks; correspondingly, the 

Q . meson, assumed to be a pion, is described as a quark-antiquark pair. Explicit expres- 

'- ^ ', sions for the unpolarized generalized parton distributions are obtained and evaluated 

Q^' in different kinematics. 
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: I. INTRODUCTION 

The fundamental role of a nonperturbative pion cloud surrounding the nucleon is well 
explained in Quantum Chromodynamics (QCD) as a consequence of the spontaneously- 
broken chiral symmetry. The pion cloud associated with chiral-symmetry breaking was first 
discussed in the context of deep-inelastic scattering (DIS) by Feynman |]J and Sullivan 
As realized by Thomas it can give an explanation of the flavor-symmetry violation in 
the sea-quark distributions of the nucleon thus naturally accounting for the excess of d 
(anti)quar 
rule 

entirely to its virtual meson cloud |8], the role of mesons in understanding these data has 
been extensively discussed in connection with parton distributions (for reviews, see Refs. 0, 



:ks over u (anti) quarks as observed through the violation of the Gottfried sum 
y, ^, 0|- Although the nuclej^n's nonperturbative antiquark sea cannot be ascribed 



IG 



ul)- 
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A similar important role is expected to be played by the meson cloud in the case of 
Generalized Parton Distributions (GPDs) that have recently been introduced and discussed 
in connection with Deeply Virtual Conipton Scattering (DVCS) and hard exclusive meson 
production (for reviews, see Refs. 



15|). Due to their intrinsic nonperturbative 



nature GPDs cannot be calculated directly within QCD. Up to now, only first results of their 
Mellin moments have been obtained by lattice calculations (see, e.g., Ref. Q]). Therefore, 
in order to guide the planning of possible experiments one has still to rely on models. 

The first calculation performed in the MIT bag model did not consider the pion 
cloud explicitly. Further calculations have been performed in the chiral quark-soliton model 



lia 
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21[ . The model is based on an effective relativistic quantum field theory where the 
instanton fluctuations of the gluon field are simulated by a pion field binding the constituent 
quarks inside the nucleon. The model is theoretically justified in the limit of the large number 
of colors Nc, and in the leading order in the l/N^ expansion it is not possible to obtain results 
for separate flavors, only special flavour combinations being nonzero. GPDs have also been 



22j |. A first attempt to explicitly include 



calculated within the Nambu-Jona-Lasinio model 

the meson cloud in a model calculation of GPDs has been discussed in Ref. [23] in terms of 
double distributions [2^ . 

A complete and exact overlap representation of GPDs has recently been worked out 
within the framework of light-cone quantization 0,12^. A Fock-state decomposition of the 
hadronic state is performed in terms of A^-parton Fock states with coefficients representing 
the momentum light-cone wavefunction (LCWF) of the A^ partons. The same approach has 



been followed in Refs. 
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where GPDs both in the chiral-even and chiral-odd sector 
were derived assuming that at the low-energy scale valence quarks can be interpreted as the 
constituent quarks treated in constituent quark models (CQMs). This assumption is based 
on the idea that there exists a scale Qq where the short range (perturbative) part of the 
interaction is negligible and, therefore, the glue and sea are suppressed, while a long range 
(confini ng) part of the interaction produces a proton composed by (three) valence quarks, 
mainly j^Q]. Jaffe and Ross 0] proposed to ascribe the quark model calculations of matrix 
elements to that hadronic scale Qq. In this way, quark models summarizing a great deal 
of hadronic properties may substitute for low-energy parametrizations, while evolution to 
larger momentum is dictated by perturbative QCD. 

In this paper, we will study the possibility of integrating meson-cloud effects into the 



valence-quark contribution to GPDs. Along the lines originally proposed in Refs. jS^, |3c 
a meson-baryon Fock-state expansion is used to construct the state |A^) of the physical 
nucleon. In the one-meson approximation the state \N) is pictured as being part of the 
time a bare nucleon, |A^), and part of the time a baryon-meson system, \BM). In this 
framework, it will be shown how to apply the convolution approach used for the standard 
parton distributions in deep inelastic scattering 3^ to the case of GPDs. The main idea of 
the convolution approach is that there are no interactions among the particles in a multi- 
particle Fock state during the interaction with the hard photon. Therefore the external 
probe can scatter either on the bare nucleon, |A^), or on one of the constituents of the 
higher Fock states, \BM). In the socalled DGLAP region, with one (anti)quark emitted 
and reabsorbed by the physical nucleon, the GPDs are obtained by folding the quark GPDs 
within bare constituents (nucleons, pions, deltas, etc.) with the probability amplitudes 
describing the distribution of these constituents in the dressed initial and final nucleon. In 
contrast, in the socalled ERBL region where a qq pair is emitted from the initial nucleon, 
the GPDs are obtained from the overlap between wavefunctions of Fock states with different 
parton number, which corresponds, in the meson cloud model, to the contribution of the 
interference terms between the \BM) component in the initial state and the bare nucleon 
in the final state. 

The model is revisited in Section ITTl where all the necessary ingredients that one can find 
scattered in the literature are consistently rederived to construct the LCWF in the meson 
cloud model. The definition of the unpolarized GPDs are introduced in Section IIIIl and 
the convolution formulas for the GPDs in the meson cloud model are explicitly derived 
in Section IIVI In Section El we describe the model calculation considering the case of a 
pure pion cloud, and assuming a constituent quark model to construct the LCWFs of the 
bare nucleon and the constituents of the baryon-meson Fock state. Numerical results are 
discussed in Section IVII and some conclusions are drawn in the final Section. Derivation of 
auxiliary quantities is detailed in three appendices. 



II. THE MESON-CLOUD MODEL FOR THE NUCLEON 



The basic assumption throughout this paper is that the physical nucleon is made of 
a bare nucleon A^ dressed by the surrounding meson cloud so that the state of the physical 
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nucleon is decomposed according to the meson-baryon Fock-state expansion as a superpo- 
sition of a bare nucleon state and states containing virtual mesons associated with recoiling 
baryons. This state, with four-momentum = {Pm,Pm,Pn±) = {Pn^Pn) and helicity A, is 
an eigenstate of the light-cone Hamiltonian 

Hlc = E [Hiiq) + <{q) + Hj{N, BM)] , (1) 

B,M 

i.e. 

r»2 + /If 2 

HlcIPn, A; N) = X \PN, A; N). (2) 

Pn 

Here H^lq) stands for the effective-QCD Hamiltonian which governs the constituent-quark 
dynamics, and leads to the confinement of three quarks in a baryon state; analogously, 
H^{q) describes the quark interaction in a meson state. Thus we assume that the three- 
and two-quark states with the quantum numbers of a baryon and a meson are the eigenstates 
ofH^iq) 3.ndH^{q),e.g. 

Hi\ps,X;B) = ^^±^\ps,X;B), (3) 
Pb 

<|PM, A; M) = ^^^^^\pM, A; M). (4) 

Pm 

In Eq. (^, Hi{N, BM) is the nucleon-baryon-meson interaction, and the sum is over all the 
possible baryon and meson configurations in which the nucleon can virtually fluctuate. Using 
perturbation theory, we can expand the nucleon wavefunction in terms of the eigenstates of 
the bare Hamiltonian Hq = H^^q) + H^^^q), i.e. 

\PnA]N) = VZ {\pn,X]N) + }^ — 

V- / \n2){n2\Hi\ni){ni\Hj\pN,\]N) \ 

where J2' indicates the summation over BM intermediate states, and Z is the wavefunction 
renormalization constant. In the one-meson approximation, we truncate the series expansion 
of Eq. (0) to the first order in Hj, and as a result we obtain 

\PN,X;N) = Vz\pr,,X;N) + J2 \pnA\N{BM)) 

B,M 

dPfid^Psx /"dp^d^pMx 
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In Eq. the normalization factor \fZ only affects the bare core |iV), and not the meson- 
baryon component. As discussed in details in Refs. ^ Is^, this prescription is consistent 
with assuming that the nucleon-baryon-meson coupling constant in if/ is taken equal to the 
renormalized value gNBM-, related to lowest order to the bare coupling g%BM 9nbm = 

NBM- 

Finally, the hadron states in Eq. © are normalized as 

{p'\ p1, A'; H\p\ p^A; H) = 2{27rfp^5{p'^ - p'')5^'\v'i. - vAw- (7) 

A. The nucleon wavefunction 

In this Section starting from Eq. (jH)) we derive the explicit general expression of the 
nucleon wavefunction on the basis of bare-nucleon and baryon-meson Fock states. 

We first evaluate the energy denominator in Eq. ® using the following expression for 
the energy of the particles in terms of light-front variables 

If we write the momenta of the baryon, pB-, and the meson, pMi in terms of the intrinsic 
(nucleon rest-frame) variables, i.e. 

vb = yPN, = (1 - y)PN, 

Pb± = k± + y pn±, Pm± = -k± + (!-?/) Pn±, 

we have 



(9) 



1/2 M2+k2 M2 +k2 

{En -Eb- Em) = A4 ^ ^ ^ 



^pn \ y ^-y ) 

^(M^-M|^,(2/,k^)), (10) 



2Piv 
where 

MsMVy^ k^) = + ^ _^ , (11) 

is the invariant mass of the baryon-meson fluctuation. 

Furthermore, the transition amplitude {B{pb, X')M{pm, A")|ii/|A^(pAr, A)) in Eq. (0) can 
be rewritten as 

{B{pB, X')M{pM, X")\Hj\N{pN, A)) = 

= i27rf5ipt+p+-pt,)5^'\pB±+PM± - PN±)V,),.iN,BM), (12) 
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where the vertex function Vy y, {N, BM) has the following general expression 

V^)^y,{N, BM) = UNaipN, A)t;"^^x/3(PM, X")iJ^{pB, A'). (13) 

Here is the nucleon spinor, x ^iid ip are the field operators of the intermediate meson 
and baryon, respectively, and a, P, 7 are bi-spinor and/or vector indices depending on the 
representation used for particles of given type. 
Using the results of Eqs. (fTUj) and (O, we find 

\P,, A; N) = ^\p,, A; AT) + E / E <P'Jr''\y, k.) 

X \yp+, + ypN±, A'; B) \{1 - y)p+, -k^ + (1 - y)pN±, A"; M), (14) 

where we introduced the function (j)'^!'y,'^^^\y,'k±) to define the probability amplitude for 
a nucleon with helicity A to fluctuate into a virtual BM system with the baryon having 
helicity A', longitudinal momentum fraction y and transverse momentum k^, and the meson 
having helicity A", longitudinal momentum fraction 1 — y and transverse momentum — k^, 
i.e. 

We note that Eq. ()14|) is equivalent to the expression of the nucleon wavefunction obtained 
in the framework of "old-fashioned" time-ordered perturbation theory in the infinite mo- 
mentum frame (see Ref. 8^). 

By imposing the normalization of the nucleon state as in Eq. (j?)), from Eq. ()14|) we obtain 
the following condition on the normalization factor Z 

l = Z + Pbm/n, (16) 

with 

_ ^ dyd^k^ 1 \Vi/UN,BM)\^ 
^^/^ 1^,7 2(27r)3 y(l-y)^f3^,,[M^-Ml^,(y,k^)p- ^'^ 

Here Pbm/n is the probability of fluctuation of the nucleon in a baryon-meson state, and, 

accordingly, Z gives the probability to find the bare nucleon in the physical nucleon. 



B. Partonic content of the nucleon wavefunction 



In this Section we derive the expression of the nucleon wavefunction ()14|1 in terms of the 
constituent partons of the nucleon core and of the meson-baryon components. 
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The light-front state of the bare nucleon is given by 



\p^,x;n) = y: J 



f 


dx 







4 = 1 



where '$x '^^\{xi, k_Lj; Aj, Ti}j=i^2,3) is the momentum hght-cone wavefunction which gives the 
probabihty amphtude for finding in the nucleon three quarks with momenta (xjp^, p^^ = 
ki± + XiPN±), and spin and isospin variables Aj and r^, respectively. In Eq. flTH|) and in the 
following formulas, the integration measures are defined by 



dx 



N 



v.=i2(2vr)3^ 

Next we consider the component of the meson-baryon Fock state in Eq. fll4p . 
The light-front state of the baryon is given by 

3 



(19) 
(20) 



\PB, A'; B)=Y.j 



/ 


dx 







[d^k_L]3^f;'^^({xi,k_Li; Ai,rj}i=i,2,3) n l^i^B' "Pi^AhThq), (21) 



J 3 



where now the intrinsic variables of the quarks x-i and kf refer to the baryon rest frame, i.e. 
= Pi IPb and pi_L = kn. + XiPsi. {i = 1, 2, 3). 

An analogous expression holds for the light-front state of the meson, i.e. 

da;4dx5 dk4_Ldk5_L 



\PM, A"; M) = Y.J 



X 



5(1 - X4 - Xs) 5^'nk4± + k 



5± 



'^y;^^\{xi,k±i;Xi,Ti}i=4^5)Y[ \xiptj, Pi±,Xi,Ti;q), 



(22) 



i=4 



with Xi = pt/pii, Pi± = XiPM± + ki_L (i = 4, 5). 

When we insert the expressions of the baryon and meson states in Eq. (fT^ . it is convenient 
to rewrite the kinematical variables of the partons as follows. 



For i = 1, 2, 3: 



Xi 



pb pnPb 

XiPB± + ki_L = Xi (k_L + y Pn±) + k^. 

^iPN± + ki_L + Xi k_L = ^iPN± + k-_L, 
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where = pf /pj^ is the fraction of the longitudinal momentum of the nucleon carried by 
the quarks in the baryon, and k'^j^ is the intrinsic transverse momentum of the quarks with 
respect to the nucleon rest frame. 
For i = 4,5: 

pi _ Pt Pn _ 6 

Xi 



Pm Pn Pm 1 - y 
Pi± = XiPM± + kj_L = Xi (-k_L + {1 -y) pn±) + 

= iiPN± + ki_L -Xikj_= ^iPN± + k-_L, 

with = pf /pij and k'-_|_ the intrinsic variables of the quarks in the meson with respect to 
the nucleon rest frame. Accordingly we transform the variables of integration as follows. 
For i = 1, 2, 3: 

Xi ^ y Xi, 

kj± k-_L = ki_L +Xik_L. 

For i = 4,5: 

Xi ^ ^i = {1- y)xi, 

k't± ki_L = kj_L Xi kj_ . 

The meson-baryon component of the nucleon wavefunction in Eq. p4|) can then be written 



as 



X 



A', A" \,n ^^^N - J-^hMKUi ^-LJ 

X 



^5f;[^l({e„k^^;A„r,},=4,5)n I^^P^^' k'^+6Piv±, A„ r,;g), (23) 



where the wavefunctions \E'f,' and \E'^f,' incorporate the Jacobian J of the transformation 
Xi ^ ^i, i.e. 



l^^f;[^l({x,,ki^;A„r,},=i,2,3), (24) 
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(1-2/) 

Finally, by introducing the following definition 



^^^'''lf\{x,, A„ r,},=4,5). (25) 



^f3;'„M2,-M2,^(l/,k^) 
X ^f; k:^, A„ r,},=i,2,3)^?; '^'({e^, k^^, A„ r,},=4,5), (26) 

Eq. can be simplified to the following expression: 

,.,A;iv(.M)) ^ /d.d^k. /;n f r g § /n ^ 

X 5 - E 6) (k^ - E k:^) 5 (i - E e.) 5^^) (e k'.^) 

X E^A''^^'(2/,k^; {e,,k:^;A„r,},=i„„,5) 

5 

X n l^iPN^ k'^ + 6p7v±, \,ri;q), (27) 

i=l 

where ^'^^ can be interpreted as the probability amplitude for finding in the nucleon a 
configuration of five partons composed by two clusters of three and two quarks, with total 
momentum {yptr,PB±) and ((1 - y)p%,PM±), respectively. 



III. THE UNPOLARIZED GENERALIZED PARTON DISTRIBUTIONS 



In the definition of GPDs it is useful to choose a symmetric frame of reference where the 
virtual photon momentum and the average nucleon momentum = \{p% + p'n) 
collinear along the z axis and in opposite directions, i.e. 



Pn 



p'n 



2(1 + ' 
2(1 -Op^ ' 



;1 + 0P7^ , 

;i-0p^, 



2 
A 



+ ' 



Pn 



+ A5_/4 

2(1 + 0^ ^ 

"m^ + A^ /4 . ' 

Pn 



2(1 -Op^ ' 



(28) 



Furthermore, = —q^q^ is the space-like virtuality that defines the scale of the process. 



t = A2 



{p'^ — p%y is the invariant transferred momentum square, and the skewness ^ 



describes the longitudinal change of the nucleon momentum, 2,^ = — A+/p 



AT- 
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According to Ref. [Ij], for each flavor q tlie soft amplitude corresponding to unpolarized 
GPDs reads 

(29) 

where x defines the fraction of the quark light-cone momentum {k^ = xp^), Ajv {^'n) 
is the helicity of the initial (final) nucleon, and the quark-quark correlation function is 
integrated along the light-cone distance z~ at equal light-cone time {y^ = 0) and zero 
transverse separation (z_l = 0) between the quarks. The leading twist (twist-two) part of 
this amplitude can be parametrized as 

1 

1 _ ^ za+^Aj, 



where H'^{x,C,, and E'^{x,^, A.±) are the chiral-even helicity conserving and helicity 
flipping GPDs for partons of flavor q, respectively. Taking different proton-helicity combi- 
nations, we have 

F|+(x,e,A^) = F^_(x,e,A^) 

= H^ix, A^) - E'^ix, e, A^), (31) 

F^(x,e,A^) = -{FUx,^,A^)r = r^y!^-^^E^{x,^,A^), (32) 



where 

+ 2A2 



(33) 



and 



2 /\/r2 

N 



- to = (34) 



is the minimal value for — t at given ^. 

IV. THE CONVOLUTION MODEL FOR THE UNPOLARIZED GPDS 



Before deriving the convolution formulas for the GPDs in the meson cloud model, we need 
to specify our conventions for the kinematical variables. The Fock expansion of the nucleon 
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wavefunction in Eq. (jl8j) as well as the LCWF for the hadron states in Eqs. (j2ip . (j22p and 
fj27|) do not depend on the momentum of the hadron, but only on the momentum coordinates 
of the constituents relative to the hadron momentum. This fact reflects the well known result 
in light-front dynamics that the centre of mass motion can be separated from the relative 
motion of the constituents. On the other hand, the arguments of the LCWF can most 
easily identified in reference frames where the hadron has zero transverse momentum. We 
call such frames "hadron frames" , and more specifically we introduce the names "hadron-in" 



and "hadron-out" for the 
momentum, respectively 



rames where the incoming and outgoing hadron has zero transverse 



25l |. In general, we denote the momenta of constituents belonging 



to the incoming hadrons with unprimed, and the momenta of constituents belonging to the 
outgoing hadron with primed variables. Furthermore, we label quantities in the hadron- 
in (hadron-out) frame with an additional tilde (hat). The relations between the momenta 
of the constituents in the "average frame" defined in Eq. (j28|) and the variables in the 
hadron frames are obtained via a transverse boost, i.e. a transformation that leaves the plus 
component of any four-vector z unchanged. It reads 



(35) 



where the two parameters and b_|_ are given by 6+ = (1 + OPiv ^^'^ t)_L = — A_|_/2 for 
the transformation from the average frame to the hadron-in frame. Likewise, a transverse 
boost with parameters = (1 — ^)p^ and b_L = +A±/2 leads from the average frame to 
the hadron-out frame. 

In the following we will derive the convolution formulas for the GPDs in the three different 
regions corresponding to < x < 1, — ^ ^ x < C,, and — 1 < x < 



A. The region ^ < x < 1 

In this region the GPDs describe the emission of a quark from the nucleon with momentum 
fraction and its reabsorption with momentum fraction x — ^. In the meson-cloud model, 
the virtual photon can hit either the bare nucleon or one of the higher Fock states. As a 
consequence, the DVCS amplitude can be written as the sum of two contributions 

^AlA.(^,e, A^) = ZF,t;;(a:,e, A^) + . (x,e, A^), (36) 
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where F'^''"^^*^ is the contribution from the bare proton, described in terms of Fock states 
with three valence quarks, and 5F'^ is the contribution from the BM Fock components of 
the nucleon state, corresponding to five-parton configurations. This last term can further be 
split into two contributions, with the active quark belonging either to the baryon [5F'^/^^^) 
or to the active quark inside the meson (6F'^^^^^), i.e. 




Figure 1: Deeply virtual Compton scattering from the bare nucleon. 

The valence-quark contribution corresponding to the diagram of Fig. ^ can be calculated 
in terms of the light-front overlap representation derived in Ref. 2^, and applied here to 
the case of = 3 valence quarks. It reads 

nCZ(^^ A^) = E E Ss,, [ [dx]3[d^kj3 Six - -X,) 

X vI/;^[^l*({x'„k',;A„r.})vI/;^[^l({5.k.;A„r.}), (38) 

where the LCWF ^^'''^'({xi, k^^; A^, Ti}) is the bare-nucleon LCWF of Eq. jlHl), and Sj 
labels the quantum numbers of the jth active parton. The integration in Eq. ()38|1 is over 
the average quark transverse momenta kj^j and the average quark longitudinal momentum 
fractions Xi = kf /p^. The kinematical variables appearing as arguments in the LCWFs 
in Eq. (jHHj) have been defined in the "hadron-in" and "hadron-out" frames following the 
conventions of Ref. j25|. In particular, the momenta of the partons belonging to the incoming 



hadron are given by 



ii = . k^i = k±i + . (39) 



for the spectator quarks {i j] 
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for the active quark. Likewise, the LCWF arguments for the outgoing hadron are exphcitly 
given by 



1-e 

for the spectator quarks {i ^ j); 

-J _ ~ ^ 

for the active quark. 



Xi 



1 2 



1 - Xj Aj 



(42) 




(a) 




Figure 2: Deeply virtual Compton scattering from the virtual (a) baryon and (b) meson components 
of a dressed nucleon. 

We now consider the contribution from the BM component of the proton wavefunction, 
starting from the SFy^^^ term in Eq. (jHTj) represented in Fig. 2a. In this case, the baryon is 
taken out from the initial proton with a fraction ys + C, of the average plus- momentum p^, 
and after the interaction with the initial and final photons is reinserted back into the final 
proton with a fraction ys — ^ of the average plus-momentum p^. The transverse momentum 
of the baryon is pb± — A_l/2 before, and p^_|_ + A_l/2 after the scattering process. The 
meson substate is a spectator during the whole scattering process, with momentum in the 
initial and final state equal to (pjj, — Pb±)- 

The baryon contribution to the 6Fy scattering amplitude can be easily evaluated by 
calculating the matrix element in Eq. ()29j) between the BM components of the initial and 
final nucleon expressed in the hadron-in and hadron-out frame, respectively. For the initial 
state, they are given by 



dye dyM ~ - ^ f d^k^^d^k 



/ -j=^=6{l -ys- yhi) / TTTT-T^ 5 (ks^ + k 



ys Vvm J 2(2 



TT) 

X E <Py>}''''''''\yB,ks^)\{yB + Ol^N,PB±-^ (43) 

A', A" ^ 
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The coordinates of the baryon and meson states in the average frame are related to the 
variables in the hadron-in frame through the transverse boost in Eq. as explained before. 
They are explicitly given by 

Vm = ' = Pm± + ^ ■ (45) 



Analogously, the BM components of the final nucleon state are given by 

2(27r) 



BM-' \ yB\ yM 



X E ^fx'^'''^'\yBMBA) K^B -Op^,Pbx + ^,A')|pl-„-p^^, A"), (46) 

A', A" ^ 

where the coordinates in the hadron-out frame and in the average frame are related by 

i/M — \ — ^ ' ~ ~ \ — ^~2r' 

The final result for the baryon contribution to the 6Fy scattering amplitude is 

^i^A,A.(-,e,A,) - ^^E^E J, ^i^VA (-'-'^-j 

X <^J;^i"'""^(y., k.,) [0&r^""^(^^, k..)]*, (49) 



where 



[Vb Vb J 2Jl - e/vl ^ 



X (PB^PB±,A'|V^(-|)7+^(^)bJ,PB±,A) (50) 

is the scattering amplitude from the active baryon in the BM component of the nucleon. 
The overlap representation of -F^a'^ terms of LCWFs is given explicitly in App. 1X1 

Analogously, we can derive the meson contribution to the scattering amplitude, corre- 
sponding to the case when the pion takes part to the interaction process while the baryon 
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remains as a spectator (see Fig. 2b). In such a case the role of the meson and baryon 
substates is interchanged with respect to the situation described before. The meson is taken 
out from the initial proton with a fraction yj^,[ + ^ of the average plus- momentum p^, and 
after the interaction with the initial and final photons is reinserted back into the final proton 
with a fraction — ,^ of the average plus- momentum p^. The transverse momentum of 
the meson is pjv/j, — A_l/2 before, and Pa/_i_ + A_l/2 after the scattering process. Viceversa, 
the baryon substate is inert during the whole scattering process, with the same momentum 
(p^, —'Pml) the initial and final state. Therefore the meson contribution to the 5Fy 
scattering amphtude is given by 

^ ,Ajv (Af,BM)/-| ~ f X r ,A' (Af,BA/). ^ .-.^ , . 

X (t^\"\ (1 - Vm, -kAf±) [0A"A' (1 - Vm^ -kAf±)J , (51) 



where 



Vm Vm ) l^l-eitu^ 

X (p'A^,p'M±,A'R-^)7+^(^)bl^,PM±,A) (52) 

is the scattering amplitude from the active meson in the BM component of the nucleon. In 
Eq. (jKT|) . the meson coordinates in the hadron-in and hadron-out frames are related to the 
variables in the average frame by 

~ _ Vm + i r _- 1 - Vm A± 

~ ~Y+Y ' ~ 1+^^ — 2~' 

_ Vm - e ^ _ I . . 

Vm — i_ ^ ' Af± — Pm± + _ ^ 2 ■ ^ 

The explicit expression of F^l^^ in terms of LCWFs is given in App. IHl 



Forward limit 



In the limit 0, where x goes over to the parton momentum fraction x, and ^ = 

A± = 0, the scattering amplitude without nucleon helicity fiip reduces to the ordinary 
parton distribution, i.e. 

F|+(x,0,0) = Fl_{x,0,0) = q{x) = Z q'"'''^{x) + 6q{x). (55) 
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In Ref. 



27|, we derived the forward limit of the valence-quark contribution, i.e. 



pq,bare^—^ 0, 0) = q^"'^^{x). Here we show that the forward limit of Eqs. and gives the 
5q{x) contribution to the parton distribution considered in Ref. |9| within the meson-cloud 
model. 

In the case of the active baryon, we have Vb = v'b = Vb and k^x = k'^x = ^b±, and 
Eq. reduces to the following expression 

^F^r (-,0,0) = E Tdy. / E \<Ptr''\yB,i^B.)\' Fli- f^,o,o) 

^Jx J 2{2'KYyB XVb j 

= H\ — ) fBM,N{yB), (56) 



M 

where the splitting function 



^' Vb Wb, 



d^k 



ibmAvb) = / ^ E \<Pt;:^r'''\yB,i^B±)\' (57) 

coincides with the definition given in Ref. P]. 

Analogously, in the case of the active meson we have yu = v'm — Vm and k^x = ^'m± — 
kM± and Eq. reduces to the following expression 



5F|r(x,o,o) = EEr^/ 



d^k 



B A,A" 



B 

where we use the definition 



Ix VM J 2(27r)3 

™a-^M,-kMx)pFTr f-,o,o' 



E / -fMB,N{yM)qM — , (58) 



JmbAVb) = [ tf^^ii E I0AA'^'^*^^(1 - VM, -^M±)f- (59) 
^ ^^^^) X,X" 

Since the probability for the dressed nucleon to consist of a bare baryon and meson is 
independent of which one interacts with the probe, we assumed in Eq. (j^^ the following 
condition 

fMB,N{yM) = fBM,N{^ - Vm)- (60) 

This property cannot be derived a priori from the definition of the convolution model, but 
it is an additional physical input which may be used to restrict the vertex functions of the 
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model as we will discuss in SectjV] Moreover, the relation automatically ensures both 
momentum and baryon number sum rules. 

As a final result, the contribution of the higher Fock states to the parton distribution is 
given by 



Sq{x) = J2 



MB L 



(61) 



which coincides with the formulation of the meson cloud model in deep inelastic process 

n 

(see, for example, Ref. [2| and reference therein). 



B. The region — 1 < x < — ^ 

In this region, the scattering amplitude describes the emission of an antiquark from the 
nucleon with momentum fraction —{x + ^) and its reabsorption with momentum fraction 
— (x — ^). As a consequence, the only non vanishing contribution can come from the active 
antiquark in the meson substate of the BM Fock component of the nucleon wavefunction, 
i.e. 

^aV.(^'^'^^) = (62) 

where SFy^^{x,^, A±) corresponds to the meson scattering amplitude illustrated in Fig. 
2(b), and is given by 

X <p',;:^''-''">(l - VM, -k«,x) - v'm, -k«,x)|*. (63) 

Here the relations between the coordinates in the hadron frames and in the average frame 
are the same as in Eqs. and (jSH). In addition, we note that Eq. corresponds to the 
same convolution formula (jSH), with the integration range over between —x and 1, and 
with the explicit LCWF overlap representation of -F^/f^ in the range — 1 < x < — ^ given in 
App.E 



Forward limit 
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The scattering amplitude from the meson substate has the following forward limit: 

Fff 0, 0) = g(^) = ). (64) 

Vm Vm Vm 

As a consequence, in the forward limit Eq. (jU^ reduces to 

5F:f ^(x,0,0) = -E /' (-—] fMsAVM) 

B J-^ Vm \ Vm J 

= -6q{-x). (65) 

C. The region — ^ < x < ^ 

In this region, the scattering amplitude describes the emission of a quark- ant iquark pair 
from the initial proton. As discussed in Ref. Q], in the Fock-state decomposition of the 
initial and final nucleons we have to consider only terms where the initial state has the same 
parton content as the final state plus an additional quark-antiquark pair. In the present 
meson-cloud model, the initial state is given by the five-parton component in Eq. ()27j) . 
while the final state is described by the three- valence quark configuration given in Eq. (fTH|) . 
multiplied by the normalization factor ^/Z . To label the coordinates of the initial and final 
partons we follow the conventions explained in Ref. 2^. This means that we use the same 
numbering for the spectator partons in the LCWFs of the initial and final state partons. 
Thus the three partons in the outgoing nucleon are numbered not as z = 1,2,3, but as 
i = 1, ...,5 with the labels of the active quark (j) and the active antiquark (j') omitted. In 
principle, we can distinguish between two cases: i) the active quark belongs to the baryon 
substate (j = 1,2,3), and the active antiquark is in the parton configuration of the meson 
substate of the initial nucleon; ii) both the active quark and antiquark belong to the meson 
substate of the initial nucleon [j = 4 and j' = 5) and the baryon is a spectator during 
the scattering process. However, this last contribution is vanishing because it involves the 
overlap of two orthogonal states, i.e. the wave functions of the baryon in the initial state and 
of the bare nucleon in the final state. As a consequence, the only non vanishing contribution 
corresponds to the case i) which is pictured in Fig. El 

The final result for the scattering amplitude in the ERBL region is derived using the 
overlap representation of Ref. j25] with the initial wavefunction describing the initial 5- 
parton configuration given in Eq. (|27|) and the final state described by the valence-quark 
wavefunction in Eq. (fTS|) . It reads 
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Figure 3: Deeply virtual Compton scattering in the ERBL region. 



1 p^i^.c^2 ^s„s,Sx„ -A,., / dye / dyM^il -ye- Vm) 

1 - U (l + U B,M A„r,i=l j' ' J J 



X / d'kB± / d'kM±5ikB± + k 



M±) 



dxi 



X 



n dx, / d^k^a / n d'k 



1=1 



[2(27r; 
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V 



^¥=3,3 



i = l 



X ^f/lfe, kB±; {i„ k,^, A„ r,},=i,...,5)[^5j^l({^^:, A., r,},=i,2,3)]*, (66) 



where is a color factor which comes from the color component of the initial- and final- 
state wavefunction. By taking the color component equal to '}2ijk^/V^Sijk\q^q'q^) and 
Y^ij ^1 ^/^ ^ijWo.'^ ) for the baryon and meson state, respectively, we have K^, = l/VS. 

In Eq. (jnni) the coordinates of the initial partons in the hadron-in frame are related to 
the parton momenta in the average frame by 

Xi A I 



Xi 



Xj-i 



k±i + 



k±i = k±j 



k±/ - -k_Lj - 



1 + e 2 ' 

1 - Xj A_L 

Try ~ ' 

1 + x^ Aj 



for % ^ 



(67) 



where for the average momentum of the spectators partons we used the standard definition 
= I {h + A;-), while for the active quark we defined kj = \ {kj — k'j), i.e. half the relative 
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momentum between the active quark and antiquark. The corresponding relations for the 
coordinates of the final partons in the hadron-out and average frames are 

^'i = ' ^'±i = - ' for i ^ . (68) 



1 - e 1 - e 2 

Finally, the five-parton component of the initial state wavefunction ^^x^^^ -^1- ^^^P is 
given by 



1 1 _ V^,\„iN,BM) 



E 



X ^f; ^^\{Ci, Ki±, Xi, ri}i=i,2,3)^5" ^^\{Ci, ^^±, ri}»=4,5) (69) 

with the coordinates in the baryon-in [Q, with i = 1,2,3) and meson-in [Q, with i = 4,5) 
frames defined as 

~ Xi ~ ~ ~ 

Ci = — , Ki± = ki_L — Ci\s.B±, ^ = 1,2, 3, 

~ % ~ ~~ (70) 

Ci = -r^ , «i± = kj_L - OkM±, i = 4, 5. 
Vm 

V. MODEL FOR THE PION CLOUD OF THE PROTON 

In this section we specify the ingredients for the model calculation of the unpolarized 
GPDs of the proton. We restrict ourselves to consider only the pion-cloud contribution 
disregarding the contributions from mesons of higher masses which are suppressed. As a 
consequence, the accompanying baryon in the l-Bvr) component of the dressed proton is a 
nucleon or a A. 

Vertex functions 

The vertex functions for the transition p Bit are given by 



^A^o(P'^7^) = i gpN^Ux'{pN)lbUx{Pp), 

^A^o(P^^^) = '^9pAnU'^'{PA){Pp-PA)l^Ux{Pp) (71) 
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where u\{Pn) and u^,{p/\) are the nucleon spinor and the Rarita-Schwinger spinor, respec- 
tively, defined in Appendix lO 

Because of the extended nature of the vertices one has to replace the coupling constants 
in Eq. (frT|) with phenomenological vertex form factors, Gatba/ (y, k\), which parametrize the 
unknown dynamics at the vertices. We use the following parametrization j^o] 

where A^m is the cut-off parameter and the invariant mass of the baryon-meson fluctuation 
M^^j is given in Eq. We note that the dipole parametrization of the vertex form factors 
satisfies the condition GjyBMiy, ^1) = Gnbm{^ ~y^k\) which automatically guarantees the 
property (lUUl) for the splitting functions. 

For the coupling constants we use the numerical values given in Refs. I22l|, i-e. 
5'ir7V7r/47r = 13.6 and fi-^ATr/^Ti" = 11-08 GeV-^, with qnnit = QppnO and gNAn = gpA++n-- 
The coupling of a given type of transition with different isospin components is obtained in 
terms of isospin Clesch-Gordan coefficients, i.e. gpmr+ = —V^gppwO, gpA°TT+ = —gpA+wo/V^ = 

gpA++T7- / V^- 

The violation of the Gottfried sum rule and fiavor symmetry puts constraints on the 
magnitude of the cut-off parameters. We use the values Abm = 1-0 GeV and Abm = 1-3 
GeV for the vrA^ and ttA components, respectively, because they give contributions to the 
u and d which are consistent with the requirement that the meson-cloud component of 
the sea-quark contribution cannot be larger than the measured sea quark and also with 
flavor-symmetry violation 361 1. 

n 

With the specifled parameters, in the case of the p Btt transition one has |37| . 

Pntt/p = Ppn°/p + Pnn+/p = 3Pp^0/p = 13%, 

Patt/p = Pa++tt-/p + Pa+7t0/p + Paott+/p = 2PA++7r-/p = 11%- 



LCWFs in the constituent quark model 



As explained in details in Ref. j22|, the LCWFs of the baryons and pion in Eqs. ^Zl\ and 
fl22|l ■ respectively, can be expressed in terms of the canonical wavef unctions, solutions of 
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the instant-form Hamiltonian for the A^-valence quarks of the hadron, through the following 
relation 

N 

X ^P({k.;/^.,r„/i,},=i,.,.,;v)n<A;(^c/(fc.)),(73) 

MiviA'JV i=l 

where \l/j^is the canonical wavefunction, and D^JJ^* {Rcf{ki)) are the Melosh rotations defined 



in Ref. j22|. In Eq. (ff!?jl . cui = \Jrnf + kf is the energy of the i-th. quark, and Mq = J2i is 
the free mass of the system of non-interacting quarks. 

In the case of the nucleon, we adopt the canonical wavefunction of the relativistic CQM 
model of Ref. 



I, given by a product of a space and a spin-isospin term which is SU(6) 
symmetric. For a more detailed discussion of the model see Ref. (38l |. 

The A is described as a state of isospin T = 3/2 obtained as a pure spin- flip excitation 
of the nucleon, with the corresponding spatial part of the wavefunction equal to that of the 
nucleon. 

Finally, the canonical wavefunction of the pion is taken from ref. ^ and reads 

vI/-'W(fci,4;/ii,/i2) = ^^^7^(|/^i|/i2|00)exp(-A;V(2/3^)), (74) 



with k = ki = — /c2, and P = 0.3659 GeV. The choice of the model from Ref. |39| is 
consistent with the hypercentral CQM we adopt for the nucleon and the A, since the central 
potential between the two constituent quarks in the pion is described as a linear confining 
term plus Coulomb-like interaction. As explained in more details in Ref. js^, the canonical 
expression (|7^ represents a variational solution to the mass equation. The phase of the 
pion wavefunction (fTIjl is consistent with that of the antiquark spinors of Ref. j^l • 



VI. RESULTS 



The general formalism developed in the preceding sections has been applied to calculate 
the unpolarized GPDs if and E'^ for the proton dressed by a pion cloud. The different 
contributions coming from the p BM fluctuations are derived from the basic p —>■ pn^ 
and p A++7r^ transitions. Due to the isospin relations between the different charged 
channels and the SU(6) symmetry of the spin-isospin part of the proton wavefunction, in 
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the region ^ < x < 1 the other contributions are obtained making use of the following 
relations: 

§ pu/mT+ — 2S F'^/p'^" ^ S F^/"'^'^ =25 F^^P'^° , 

^_pM/A++7r- ^ |^_pn/A+7rO^ §pu/TT°p ^ ^pd/Tx^p^ (75^ 

^pu/n+n _ /^^pu/nOp^ ^pu/w+AO _ ^pu/TT°l\+ ^ 
^pd/TT-A++ ^ 3^_pd/7rOA+_ 

Similarly, in the — 1 < x < — ^ region we have the following relations: 



(76) 



In addition, for ^ < x < 1, we have 



= (77) 

The multidimensional integration required for the numerical computation of the different 
contributions to Fy was implementec in a parallel computation using the parallelized 
version of the VEGAS routine of Ref . 41| . In this way one makes easier an otherwise time- 



consuming computation. The results presented in this Section have been obtained for some 
combinations of t and ^ as an example of the effects introduced by the sea. 
First let us study the forward limit, ^ = 0, t = 0. 

In Fig. m the spin-averaged H"^ and the helicity-flip E'^ GPDs are plotted together with 
the separated contributions from the bare proton (dashed-dotted line), the baryon (dashed 
lines) and the meson (dotted lines) in the baryon-pion fluctuation. All these contributions 



add up incoherently to give the total resu 
same as that already calculated in Ref. 



t (full curves) . The bare proton contribution is the 
27l |. rescaled by the wavefunction renormalization 
constant Z. It is always positive within its support (0 < x < 1) with the exception of E'^ 
for which it is negative. The same behaviour characterizes the baryon contribution from the 
baryon-pion fluctuation that is also limited to the range < x < 1, consistently with the 
assumption that the only active degrees of freedom for such a baryon are the valence quarks. 
Both contributions vanish at the end points of their support. The sea-quark contribution, 
extending all over the full range — 1 < x < 1, is determined by the antiquark residing in the 
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Figure 4: The different contributions to the spin-averaged {H'^, upper panels) and heUcity-flip 
lower panels) generalized parton distributions calculated in the meson-cloud model for flavours u 
(left panels) and d (right panels), at ^ = and t = 0. Dashed lines: baryon contribution from 
the \BM) component. Dotted lines: meson contribution from the \BM) component. Dashed- 
dotted lines: contribution from the bare nucleon. Full curves: total result as a sum of the different 
contributions. 

meson of the baryon-pion fluctuation. The resulting effect of the pion cloud is thus to add a 
contribution for negative x and to increase the magnitude of the GPDs for positive x with 
respect to the case of the bare proton. In particular, for positive and small x the pion cloud 
contribution as a whole is comparable to that of the bare proton, confirming the important 
role of the sea at small x found within the chiral quark-soliton model [lal20,|2l|. Since for 
X > the contribution of the baryon-meson fluctuation is effective only at low x values, the 
same faster fall off of E'^ with respect to H'^ for x — > 1 is obtained as in Ref . j^l . This is a 
consequence of the decreasing role of the Melosh transform to generate angular momentum 
in E*^ with increasing quark momentum. 

In all cases at x = the GPDs have a zero. This is due to fact that in the overlap 
integrals the various terms of the proton wavefunction are taken at one of their end points. 



This peculiar feature was discussed in Ref. 4^ and explained as an artifact due to the 



truncated Fock-state expansion. The singular behaviour of the usual parton distributions 
for X — > cannot be obtained from any finite number of Fock-state contributions, all of 
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which vanish at x = 0. 

In any case, the forward hmit of the first moment sum rules for the spin-averaged GPDs, 

rfxi/"(x,0,0) = 2, dxH'^ix, 0,0) = 1, (78) 

is correctly fulfilled. For the helicity-fiip GPDs the first moment sum rule reads 

dxE''{x,0,0) = K\ (79) 

where (k'^) is the anomalous magnetic moment of the u (d) quarks, with + k'^ = 
3(k^ + k") and — k''' = kP — /t", and being the proton and neutron anomalous 
magnetic moments, respectively. 

Experimentally, we have = 1.793 and = —1.913, thus giving = 1.673 and = 
—2.033. Without the pion cloud the values of the nucleon anomalous magnetic moments 
were found to be rather far from the experimental ones, i.e. = 0.91 and = —0.82, 
corresponding to = 1.0 and = —0.74 l27|l. This result is, however, in agreement with 

nnn 

analogous light-front calculations with point-like quarks j42l, |4^, |4q| and the values derived 



in the forward limit of GPDs derived in the Nambu-Jona-Lasinio model j22]. Including the 
pion cloud we have = 1.14 and k'^ = —1.03, corresponding to = 1.10 and = —1.07, 
closer to the phenomenological values. Similar effects are also expected for the pion-cloud 
contribution to the form factors at finite value of t, improving the agreement between the 
CQM predictions and the experimental results 

In Fig. the same results plotted in Fig. |3] are reorganized to show the isoscalar u + d 
and isovector u — d combinations. In the large- A^^c limit if" + H'^ and E'^ — E'^ appear in 
leading order in N^, while if" — H'^ and i?" + E'^ appear in subleading order Also in 
our model if" + H'^ and i?" — E'^ are found to be larger than if" — H'^ and E^ + E'^. The 
behaviour of the isoscalar combinations is also very similar to that provided by the chiral 
quark-soliton model Q, |2^ . The sea-quark contribution is an odd (even) function of x for 
the isoscalar (isovector) combination. As discussed above, the isovector combinations vanish 
at 3; = due to the truncated Fock-state expansion of the proton wavefunction. Would the 
dip at X = be filled by including the whole series expansion, the behaviour of the isovector 
combinations would resemble that derived in the chiral quark-soliton model j^, where 
GPDs in the neighbourhood of x = are mostly determined by the Dirac sea. In our model 
the sea contribution comes only from the antiquark present in the one-pion component of 
the cloud. 
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Figure 5: Isoscalar (u + d, left panels) and isovector (u — d, left panels) combinations of the spin- 
averaged (upper panels) and helicity-flip (lower panels) generalized parton distributions calculated 
in the meson-cloud model, at = and t = 0. 

Anyway, since there are no gluons in the model the total momentum of the proton is 
carried by quarks and antiquarks only. Therefore the second moment sum rule, 

J^^dxx{H'' + H'^){x,0,0) = M^ = 1. (80) 

and the spin sum rule, 

J^^ dx X (if" + H'^ + E'' + E^) (a;, 0, 0) = 2 = 1. (81) 
are consistently fulfilled in the model. In addition, 

j^^ dxx (E" + E^x, 0, 0) = 2J« - = 0. (82) 

Going beyond the forward limit we first discuss the t dependence of the GPDs at ^ = 
by comparing results in Fig. 0] at t = with those in Fig. IHlat t = —0.2 GeV^. The relative 
contribution of the different components is not modified by switching on the momentum 
transfer t, only the overall magnitude is decreased. This is in agreement with the common 
believe that the main part of the t dependence of the GPDs is exhibited by their first 
moments, i.e. by the quark Dirac and Pauli form factors. 

With a nonvanishing ^ one can explore the ERBL region with |x| < C, where in our 
model only transition amplitudes between the bare-proton and baryon-meson components 
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Figure 6: The different contributions to the spin-averaged {H'^, upper panels) and heUcity-flip 
lower panels) generalized parton distributions calculated in the meson-cloud model for flavours u 
(left panels) and d (right panels), at ^ = and t = —0.2 GeV^. Line style as in Fig. |2 

are contributing. The combined dependence of the isoscalar and isovector combinations of 
GPDs on ^ and t is shown in Figs. The t dependence at constant ^ = 0.1 can be 

extracted from Figs. Eland |Hl while the ^ dependence at constant t = —0.5 GeV^ is deduced 
from Figs. |H1 and El 

From the isoscalar and isovector combinations of GPDs plotted in Fig.[7|at ^ = 0.1 and t = 
—0.2 GeV^ we see that GPDs in the ERBL region are rather regular functions over the whole 
range, with zeros at the endpoints x = This result is quite different from the oscillatory 
behaviour predicted by the chiral quark-soliton model where the valence contribution of 
the discrete level is a smooth function extending into the ERBL region and adding to the sea 
contribution. Here this is forbidden because the support of the valence contribution is limited 
to the DGLAP region. In addition, the transition amplitude between the bare-proton and 
the baryon-meson components vanishes at the boundary of the ERBL region. As discussed 



in Ref. jl5|, the points x = ±^ correspond to very peculiar parton configurations involving 
one parton with vanishing momentum in the initial or final state hadron. As a consequence, 
when approaching e.g. ^ from below, one probes a quark-antiquark pair with one momentum 
fraction finite and the other going to zero, a configuration similar to the one of a meson 
distribution amplitude at its endpoints. On the other hand leading-twist GPDs must be 
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Figure 7: Isoscalar (u + d, left panels) and isovector (u — d, left panels) combinations of the spin- 
averaged (upper panels) and helicity-flip (lower panels) generalized parton distributions calculated 
in the meson-cloud model, at = 0.1 and t = —0.2 GeV^. Line style as in Fig. [l] 

continuous at x = ±^ to avoid logarithmically divergent scattering amplitudes in DVCS. 
In fact, this condition is fulfilled here because also approaching x = ±^ from the DGLAP 
region in our model GPDs are constrained to go to zero. This generates a discontinuity of 
the first derivative of GPDs at x = which, however, is not in contradiction with general 
principles. Including higher Fock-states of the type suggested in Ref. j42| and/or considering 
evolution to a higher scale will fill up the zero at x = 

In the DGLAP region both for x > ^ and x < — ^ no striking difference arises in Fig. U\ 
for the spin-averaged GPDs H^^"^ with respect to the results in the forward limit shown in 
Fig. El while for the helicity-flip GPDs the (negative) d contribution coming from the baryon 
in the \BM) component is responsible for a broader shape at x > ^. The t dependence of the 
different and opposite contributions to -E" and E'^ is also responsible for the small size and the 
oscillatory behaviour of at t = —0.5 GeV^ (Fig. |H1). Increasing C,, therefore compressing 
the support for the valence contribution, this effect is even more visible producing a negative 
as in Fig. IHl while the behaviour in the ERBL region remains the same. 
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Figure 8: Isoscalar (u + d, left panels) and isovector (u — d, left panels) combinations of the spin- 
averaged (upper panels) and helicity-flip (lower panels) generalized parton distributions calculated 
in the meson-cloud model, at = 0.1 and t = —0.5 GeV^. Line style as in Fig. [l] 

VII. CONCLUSIONS 

The convolution model for the physical nucleon, where the bare nucleon is dressed by its 
virtual meson cloud, has a long and successful history in explaining properties such as form 
factors and parton distributions. In this paper it has been revisited and applied for the first 
time to study GPDs. A light-front wavefunction overlap representation is obtained in the 
one-meson approximation by inserting a Fock-state expansion involving a bare nucleon and 
meson-baryon states. The model fulfills the support condition and general sum rules such 
as the number, momentum and angular momentum sum rules. Explicit expressions for the 
unpolarized GPDs have been derived and applied to the case of the meson being a pion and 
the baryon being either a nucleon or a A. 

This meson-cloud model gives the possibility to link GPDs calculated in the light-front 
formalism to the nucleon description in terms of constituent quarks including a sea con- 
tribution already at a low-energy scale and providing a suitable input for the evolution to 
higher scales. The results presented in this paper for different kinematics show an important 
contribution of the meson cloud at low x and a smooth contribution of the sea in the ERBL 
region and for negative x. As an effect of the truncated Fock-state expansion, characteristic 
nodes occurs at the endpoints of the DGLAP and ERBL region, i.e. at x = in the forward 
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Figure 9: Isoscalar (u + d, left panels) and isovector (u — d, left panels) combinations of the spin- 
averaged (upper panels) and helicity-flip (lower panels) generalized parton distributions calculated 
in the meson-cloud model, at = 0.3 and t = —0.5 GeV^. 

case and at x = ±^ in the off-forward case, where the wavefunction has to vanish. However, 
this artifact will disappear under evolution of GPDs to higher scales. In addition, since 
the contribution in the ERBL region is vanishing in the forward limit, it can not be eas- 
ily inferred from parametrizations in terms of parton distributions. Therefore, the present 
calculation gives new insights to model the off-forward features of the GPDs, and can be 
further used as input at the hadronic scale to study the behaviour under evolution at higher 
scales. 
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Appendix A: BARYON SCATTERING AMPLITUDE 



In this section we give the LCWF overlap representation of the baryon scattering am- 
phtude given in Eq. (j^UI) . and evaluated in the nucleon average frame where the plus and 
transverse components of the momentum of the initial and final baryon are given by 



This scattering amplitude gives a non-vanishing contribution only in the region ^ < x < 
1, where it describes the emission of a quark from the baryon with momentum fraction 
x + ^/y^ of the average plus momentum p^, and its reabsorption with x — ^/y^. By analogy 
with the nucleon case, we introduce the names "baryon-in" and "baryon-out" for frames 
where the incoming and outgoing baryon has zero transverse momentum. The fraction of 
plus momentum and the transverse components of the momentum of partons before the 
scattering process in the baryon-in frame will be denoted with and /Cj_L, respectively. The 
corresponding quantities for the final partons in the baryon-out frames are defined by Q and 
/Cj_L. Furthermore, defining the average momentum of the partons as ki = {ki + k'^)/2, we 
introduce the fraction C,^ = k^ /p'^ = Xi/y^. When the active parton is taken out from the 
baryon it carries a fraction + ^/y^ of ^^e average plus momentum p^, and its transverse 
momentum is kj_|_ — A_|_/2. The arguments of the initial-state LCWF in the baryon-in 
frame are obtained from the variables of the partons in the average frame by means of the 
transverse boost in Eq. (jSHI) with parameters b_L = p^^ — A±/2 and fo"*" = Psi^ + ^/Vb)- 
Furthermore, by using the spectator condition, k[ = ki = ki, one obtains 



In the final state, the active parton has a fraction of the average plus momentum of 
the baryon equal to Q — ^/ys-, ^i-nd a transverse momentum kj_L + A_l/2. The arguments 
of the final-state LCWF in the baryon-out frame are obtained from the variables of the 
partons in the average frame by means of the transverse boost in Eq. (fH^ with parameters 




(A.l) 



l+i/VB ' 

l+i/VB 




for i ^ j 



(A.2) 
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b± = Pb± + A_l/2 and 6+ = — ^/y^) and are given by 

C = 1 — ' «i± = - ^ — ( Pb + ^1 , ioii^j, 

l-i/VB ' ^'^^l-i/VB V'^^ 2 ) ■ 

Using a similar procedure as in the calculation of the valence-quark contribution to the 
scattering amplitude in the region ^ <x< 1, one finds that the LCWF overlap representa- 
tion of F^/^ is given by 

x*f/[^l*({C',K:^;A„r,,})vI/f[^]({0,K,^;A„r,}). (A.4) 



Appendix B: MESON GENERALIZED PARTON DISTRIBUTIONS 

In this section we derive the LCWF overlap representation of the meson scattering am- 
plitude in Eq. ()52|1 . with the momenta of the incoming and outgoing meson in the average 
frame given by 



Vm) 2 



Pm = Pm ( 1 - ^ ) , Pm± = Pm± + (^-^^ 

In the following we separately discuss the contributions in the DGLAP region for positive 
and negative x. 

Region ^ < x < 1 

The derivation of F'^/^^ in terms of meson LCWFs goes along the same lines as in the 
case of the baryon discussed in App. El The final result is given by 



x^Jf'[^]*({C',/c:^;A,,r,})vI/f[^l({4K,^;A,,r,}), (B.2) 



where {Q}, {K-i±} are the coordinates of the initial partons in the meson-in frame. They are 
related to the coordinates in the average frame by the same transformation as in Eq. ()A.2|) . 
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with the substitution — > yj^f, and pb± — > Pm±- Likewise, in Eq. (jR2l) {(,'}, {k[^} are 
the coordinates of the final partons in the meson-out frame and satisfy the same relations 
as in Eq. ()A.3|) with y^ and pb± replaced by and pm±, respectively. 

Region —l<x<—C, 

In this region, the meson contribution to the scattering amplitude is given by 

X vl>Jf'[/] * ({C:, k[^; A„ r,})vl/f ({0, k,r, K, r.}) (B.3) 

where the arguments of the LCWFs are given by the same expression as in the region 
^ < X < 1. 



Appendix C: VERTEX FUNCTIONS 

In this Appendix we work out the case of the — > A^vr and A^ —>■ Att transitions in 
the light-front formalism. Vertex functions for such transitions can be found in several 



places (see, e.g., Refs. 



50|), but it is convenient to show their derivation. 



The light-front vectors are defined as 



= Aj 



(C.l) 



with 

A^ = A'^±A^, A^ = {A\A^). 

We also use the notations Ar^l = A^ ± iJ^ and A = {A~^, A±). 
The light-front nucleon spinors U\{p) are given by 



Ul/2iP) 



PR 
— m 



U-l/2{p) 



72^ 



The gamma matrices are defined as in Ref. 



. 3- 



' -P. ^ 

p'^ + m 

PL 

ym — p^ J 



(C.2) 



(C.3) 
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A similar expansion for the A field involves the Rarita-Schwinger spinors given by 



W3/2(P) 


= e+i(p)^i/2(p), 




Wl/2(P) 




^e+i(p) M-i/2(p) 




= ^e(;(p)M_i/2(p) + 


^e^i(p)«i/2(p) 




= ^''-l{P)U-l/2{P), 





(C.4) 



where the polarization vectors are given by 



1 /pi -"^^ + \ 

eo P = — T ,P ,P± , 

^.(P)^(V2^,0,(-L, --!.)). (C.5) 

In the evaluation of the transition amplitudes we need the following expressions of scalar 
products 



{e+iy{pA){pN -Pa] 



kL 

1 

1 



Kr{Pi^){VN - pa), = (kl -Ml + y^Ml) , 

{e''-,nPA){PN-PA), = (C.6) 

where the relation pA± = k_L + ypN± has been used. 

The final results for the transition amplitudes Bn with nucleon helicity A = | are 

given in Table I. The corresponding results for the vertex functions with nucleon helicity 
A = — I are given by 

V^(^'^'^)(y,kx) = (-l)^/^-^Vj^^'^^^(l/,kx), (C.7) 
where kj_ = {kx, —ky). 

Our calculation for the vertex functions is in agreement with the results of Ref. js^. 



However, we found discrepancies with the results reported in Ref. 



These differences in 



the vertex functions do not affect the splitting functions of Eqs. fjHTjl and which enter 
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A' 


V{N, Ntt) 


V{N, Att) 


13 
2 2 

11 
2 2 

1 1 

2 2 

1 3 

2 2 


Mjv(l-3/) 
igNNTT ^ 


-iQNA^^iMA+yMN) 

igNA. ^ImJ^H^Ma + vMn) + (Ma + yMNf{yMN - Ma)] 
-igNAn^^^J^l + (^A + yMN){yMN - 2Ma)] 



Table I: Vertex functions for A'^ — > Ntt and N ttA. 



in the convolution formulas of parton distributions, but are important for the transverse- 
momentum dependence of the probability amplitudes (pw which enter in the convolution 
formulas for the GPDs. 
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